Interval-valued intuitionistic fuzzy multi-attribute group decision-making method considering risk preference of decision-makers and its application

An improved interval-valued intuitionistic fuzzy multi-attribute group decision-making method considering the risk preference of decision-makers is proposed to solve the multi-attribute group decision-making problem with interval-valued intuitionistic fuzzy numbers and the condition that the attribute weight information is completely unknown. Firstly, the decision-maker weight of each attribute is determined by combining similarity and proximity. In order to consider the influence of the decision-maker's risk preference on the decision result and avoid the asymptotic behavior of interval-valued intuitionistic fuzzy matrix, the risk aversion coefficient of the decision-maker is introduced and combined with the determined decision-maker's weight aggregation to form a group decision matrix. Then, the information of group decision matrix is mined, and the interval-valued intuitionistic fuzzy entropy is used to determine the attribute weight and relative weight. Based on the interval-valued intuitionistic fuzzy distance measure formula and the TODIM method, the overall superiority of each scheme relative to other schemes is obtained by calculating the superiority between schemes, and the optimal scheme is determined by comparing and sequencing. Finally, the rationality and effectiveness of the proposed method are verified by an example of mechanical assembly supplier selection decision.

and determined the weight of each attribute based on fuzzy entropy. With incomplete information on attribute weights, Wan and Dong 8 proposed a bi-objective linear programming model that satisfies both non-consistency minimization and consistency maximization to determine attribute weights. Yu et al. 9 developed a new and unified intuitionistic fuzzy multi-objective linear programming model for such portfolio selection problems. In case of completely unknown attribute weights, Jin et al. 10 defined the information entropy based on interval-valued intuitionistic fuzzy number (IVIFN) to determine the attribute weights. Xu and Shen 11 proposed an entropy measure for IVIFS and establish an entropy weight model, which is then used to determine the objective attribute weights of the alternatives. Regarding the determination of decision-maker weights, Zhang and Xu 12 proposed a goal programming method based on incomplete decision-maker weight information to determine decisionmaker weights. Wan and Dong 13 proposed a method to determine the weights of decision-makers based on similarity. The above research methods have no difference on different attributes of decision-makers' weights, which is not in line with the reality, because in practical decision-making problems, every decision-maker is often only good at certain fields but not all fields. Therefore, it is more realistic and reasonable to assign the decisionmaker weight for each attribute to each decision-maker. In summary, there is a lack of detailed and in-depth discussion on how to set the weight of decision-makers objectively and reasonably.
In the research on the similarity measure of IFS, Khalid and Abbas 14 defined a similarity measure of intuitionistic fuzzy soft sets based on Hausdorff distance, and extended this application to medical diagnosis. Zhang and Yu 15 proposed the similarity measure and the distance measure of IVIFS, which considered all the information on IVIFS and effectively overcome the information loss. Joshi and Kumar 16 studied the decision-making problem of IVIFS by calculating the similarity between each scheme and ideal scheme and extending Hamming distance. In the optimization of quantization methods for decision analysis, Behret 17 gave the conditions for additive consistency and multiplicative consistency of intuitionistic fuzzy preference relations by using the relationship between interval fuzzy sets and IFS, and established an optimization model to obtain the priority of decision schemes. Pang and Song 18 proposed a hybrid weighted aggregation method by fusing the objective comprehensive weights of experts and the weights of individual comprehensive evaluation values based on similarity, so as to obtain the group comprehensive evaluation values of schemes, and realized the ranking of schemes by defining the expected values with risk attitude factors and accurate functions. To derive priority weights of alternatives, Wan et al. 19 established an Atanassov intuitionistic fuzzy programming model, and solved by three approaches considering decision makers' different risk attitudes. Zhang 20 described the multi-stage bilateral matching problem under the intuitionistic fuzzy preference information, and constructed a stage matching weight optimization model based on the intuitionistic fuzzy similarity. Zhang et al. 21 constructed the evaluation criterion system of equipment supplier selection based on military supply chain, and proposed a group decision method based on intuitionistic fuzzy entropy and extended multi-criteria compromise solution framework.
To sum up, the above methods are very effective in solving multi-attribute group decision-making (MAGDM) problems based on IVIFN, but they also have some limitations, because different decision-makers often have different risk preferences in actual decision-making. Shi and Zhang 22 proposed an improved method of additive consistency of intuitionistic fuzzy preference relations based on decision error transfer formula for decisionmaking problems in the context of intuitionistic fuzzy preference information. Zhang et al. 23 based on the objective situation that different experts have different risk attitudes towards the same decision-making problem, introduced a risk preference coefficient to describe the different preferences of decision-makers for the uncertainty information, and constructed a ranking method from the perspectives of decision-makers' risk attitudes and similarity measures of IFS. Wang and Wan 24 investigated the group decision making with interval-valued intuitionistic fuzzy preference relations, and proposed a new order relation to rank interval-valued intuitionistic fuzzy values. Garcez et al. 25 put forward a new hybrid Grey Additive-Veto Model (GAVM) for selecting suppliers, which optimizes the choice by the decision-maker's preference. Decision-makers with different risk preferences for the same MADM problem may make different decisions, so it is necessary to consider the risk preferences of decision-makers in the process of MADM. However, the multi-attribute scheme ranking method proposed in the reference 11 does not consider the impact of the risk preferences of decision-makers on the decision results. In addition, although the fuzziness of decision-making environment is considered in the above references 17,20,21 , the psychological behavior of decision-makers is seldom considered, and some studies, although involving the psychological behavior of decision-makers, do not consider the risk attitude of decision-makers, which affects the validity and accuracy of decision-making results. In view of the fact that TOmada de decisão interativa multicritério (TODIM) method 26 is a decision-making method which is close to the preference of decision-makers on the basis of considering the psychological behavior of decision-makers, the combination of TODIM method in this study can better reflect the subjective risk preference of decision-makers, without giving the information of decision reference points in advance, and can retain the decision-making information more completely.
Global economic integration and the development of information technology also make enterprises face more intense market competition. Supply chain management, as a new management mode adapted to global manufacturing and diversified customer needs, has been widely used in enterprises. If one or several companies in the supply chain experience production blockages, changes in delivery dates, and increases in costs, the responsiveness of the entire supply chain will deteriorate, and the total cost of the supply chain will also increase.
The issue of supplier selection is an important part of supply chain management and a prerequisite for doing a good job in supply chain management 27 , which has always attracted the attention of theoretical and practical management workers, because the selection of excellent suppliers is of great significance for enterprises to reduce costs, reduce risks and enhance market competitiveness. According to the instance analysis given in this paper, the decision-making problem of mechanical assembly suppliers is essentially the research of MAGDM under uncertain conditions. Therefore, in this study, considering the fuzziness, complexity and urgency of the above-mentioned decision-making problems, a decision-making model based on IVIFS and TODIM method is established by combining the fuzziness of decision-making environment with the psychological behavior and risk www.nature.com/scientificreports/ attitude of decision-makers, and the robustness of the model is verified by different interval-valued intuitionistic fuzzy distance measurement formulas, so as to improve the scientificity and rationality of MAGDM under the risk preference state of decision-makers. Therefore, the supplier selection problem is studied as a group decisionmaking problem in interval-valued intuitionistic fuzzy environment, which is closer to the actual decision-making background, the research results of this paper have certain theoretical value and practical application value.

Preliminaries
Interval-valued intuitionistic fuzzy numbers and operational rule. Definition 1 If X is a non-empty set, then is an intuitionistic fuzzy set, in which µ A (x) and ν A (x) are respectively membership and non-membership of element , indicating that the hesitancy degree or uncertainty of element x in X belonging to A 28 .

Definition 2
If X is a non-empty set, then is an intuitionistic fuzzy set, in which μÃ(x) ⊂ [0, 1] , νÃ(x) ⊂ [0, 1],x ∈ X and Interval numbers μÃ(x) and νÃ(x) are respectively membership and non-membership of element x in X belonging to Ã .
Then IVIFS Ã can be recorded as:
) are two random IVIFN, the distance between them is expressed as 33 :

Determination of decision-makers' weights and solution of attribute weights
Determination of decision-makers' weights for attributes. Due to the differences in background knowledge and information symmetry among decision-makers, as well as the complexity and uncertainty of decision-making problems, preference conflicts inevitably exist among decision-makers. Because different decision-makers have different knowledge structures, personal preferences, and different understanding of schemes and attributes, the evaluation results will also be very different, so how to determine the appropriate weights of decision-makers has become the key to decision-making.
For each attribute C j , the evaluation values of all schemes given by the decision-maker E k about the attribute are expressed as interval-valued intuitionistic fuzzy vector r k j = (r k 1j ,r k 2j , . . . ,r k mj ) , where ω k j represents the decision-maker's weight about the attribute C j of decision-maker E k . In order to determine the value of the decision-maker's weight, two aspects need to be considered at the same time, one is similarity, which is used to indicate the similarity between the individual decision matrix of decision-maker E k and the group decision matrix composed of all decision-makers, and the other is proximity, which is used to indicate the proximity between the individual decision matrix of decision-maker E k and the group decision matrix composed of all other decision-makers except the decision-maker E k 36 . According to the above analysis, the main steps to solve the decision-maker weight ω k j of each attribute are as follows: 1. Determine the positive ideal solution vector r j under attribute C j , namely r j = (r 1j ,r 2j , . . . ,r mj ).
As we all know, most group decisions usually take the average evaluation values of multiple decision makers as the final group decision result. Therefore, the closer the evaluation value given by a decision maker is to the average value, the better the evaluation value is. On the contrary, the farther away it is from the average value, the worse the evaluation value is. Hence, this study selects the interval-valued intuitionistic fuzzy mean (IVIFM) provided by all decision makers under the attribute C j as the positive ideal solution vector r j . For the attribute C j , the evaluation value of all schemes given by the decision maker E k with regard to the attribute is expressed as the interval-valued intuitionistic fuzzy vector r k j = (r k 1j ,r k 2j , . . . ,r k mj ) , and the positive ideal solution vector r j of the attribute C j is the arithmetic mean of r k j (k ∈ K) 37 .
(3) 2. Determine all negative ideal solution vectors under attribute C j .
The negative ideal solution vectors about the attribute C j include an individual negative ideal solution vector r c j , an individual left negative ideal solution vector r L− j and an individual right negative ideal solution vector r U− j , which are expressed as: and d(r k j ,r U− j ) by using formula (4).

Calculate the similarity
The similarity between individual decision matrix of decision-maker E k and group decision matrix under attribute C j is expressed by s k j 38 : where, j = 1, 2, . . . , n ; k = 1, 2, . . . , s.

5.
Calculate the proximity ξ lk ij and average proximity η(r l j ,r k j ) between r l ij and r k ij .
6. For attribute C j , the proximity η k j between the decision-maker E k and all other decision-makers is: 7. Calculate the weight of decision-makers for each attribute. Combining the similarity and proximity obtained above, the combination weight ω k j of decision-maker E k under attribute C j is constructed.
is the control coefficient and can be weighed by changing the control coefficient between the similarity and proximity. In practical application, = 0.5 can be taken as a compromise value.
For the combination weight ω k j , the weight ω k j of the decision-maker E k under the attribute C j can be obtained after standardization: Solution of unknown attribute weights. In the process of MAGDM, attribute weight information is often completely unknown or incomplete, due to time constraints, lack of relevant knowledge and limited knowledge of the increasingly complex and uncertain group decision-making environment. In view of the fact that the attribute weight information is completely unknown in this study, the group decision matrix informa- www.nature.com/scientificreports/ tion is deeply mined, and the weight w j and relative weight w ′ j of each attribute are calculated by using intervalvalued intuitionistic fuzzy entropy.
Interval-valued intuitionistic fuzzy entropy is defined as: The relative weight of attribute w ′ j is: When the attribute information data of each scheme is an IVIFN, from the perspective of reflecting the original decision information, the more fuzzy or uncertain the attribute information is, the less information the attribute has available for the scheme is, and the larger the entropy value is, the smaller the weight should be given, and vice versa 39 . Therefore, interval-valued intuitionistic fuzzy entropy can be used to determine the weight of attributes to not only reduce the loss of evaluation information but also truly reflect the wishes of decision-makers.

Interval-valued intuitionistic fuzzy multi-attribute group decision-making method considering risk preference of decision-makers
Description of supplier selection decision-making problem. In the MAGDM problem of supplier selection, a joint evaluation team composed of production, procurement, logistics and other departments of the enterprise assigns values to the subjective and objective attributes of the selected mixed evaluation information, and ranks and selects the best suppliers according to the qualitative evaluation information given by each decision-making department and the quantitative evaluation information obtained through statistical analysis.
The evaluation value of the decision-maker E k on the attribute C j of the scheme A i is the IVIFN, which is expressed as . The intervalvalued intuitionistic fuzzy decision matrix composed of the above evaluation values is written as R k = (r k ij ) m×n .
Algorithm flow. Based on the above analysis, the flow of MAGDM method based on IVIFN is as follows: Step 1: All decision-makers E k give their corresponding interval-valued intuitionistic fuzzy matrices R k = (r k ij ) m×n .
Step 2: According to formulas (14)-(15), the decision-maker weight ω k j of the decision-maker E k about each attribute C j is calculated.
Step 3: The n-dimensional IVIFWM operator 40 is used in combination with the decision-maker weight determined in Step 2 to aggregate the single decision matrix R k = (r k ij ) m×n of the decision-maker and get the group decision matrix R = (r ij ) m×n .
It is obvious that the aggregated group decision matrix is still an interval-valued intuitionistic fuzzy matrix.
Step 4: According to formulas (17)-(18), the weight w j and relative weight w ′ j of each attribute are calculated by using interval-valued intuitionistic fuzzy entropy.
Step 5: The dominance ϕ j (A i , A k ) of scheme A i relative to scheme A k under attribute C j is calculated.
.  (20), S ij (α) and S kj (α) are the score functions of schemes A i and A k under attribute C j , respectively, and parameter θ is the risk aversion coefficient, which can be changed according to the preference of decision-makers, ranging 0 ≤ θ ≤ ( n j=1 w ′ j )/w ′ j , and the smaller the value of θ , the higher the degree of risk aversion of decision-makers 41 .
Step 6: According to the dominance matrix, the overall dominance of scheme A i relative to scheme A k is calculated.
Step 7: The overall dominance P(A i ) of scheme A i relative to other schemes is comprehensively calculated, and the schemes are sequenced according to their values to get the comprehensive evaluation values of each scheme and to select the optimal scheme.   Table 3. Interval-valued intuitionistic fuzzy matrix given by decision-maker E 3 .

Instance analysis
A regional manufacturing company is looking for the best supplier in the world for the most critical components in its assembly process. There are four suppliers A i (i = 1, 2, 3, 4) to choose from. An expert group composed of four experts (decision-makers) E k (k = 1, 2, 3, 4) from various strategic decision-making fields evaluates suppliers A i with the following five evaluation indicators (attributes): product price (C 1 ) , product quality (C 2 ) , supplier maintenance level (C 3 ) , supplier information (C 4 ) and risk factors (C 5 ) . Decision-makers E k use the IVIFN r k ij to describe the characteristics of each supplier A i under attributes C j (j = 1, 2, 3, 4, 5) , and all decision-makers E k give their corresponding interval-valued intuitionistic fuzzy matrix R k = (r k ij ) m×n . The objective evaluation information obtained through statistical investigation is shown in Tables 1, 2, 3 and 4 42 , and the optimal supplier is determined.
Step 1: Calculate the decision-maker weight ω k j for each attribute C j of decision-makers E k . For example, the decision-maker weight ω 1 1 for each attribute C j of decision-maker E 1 .
1. Calculate the similarity of decision-maker E 1 under attribute C 1 According to formula (6), the ideal solution vector r 1 under attribute C 1 is: According to formulas (7) Then, the distance between the decision matrix of the decision-maker E 1 and each ideal solution vector is obtained by using formula (4): The similarity between the decision matrix of the decision-maker E 1 and the group decision matrix composed of all the decision-makers under the attribute C 1 is calculated by using the formula (10) s 1 1 = 0.9082.

Calculate the proximity of decision-maker E 1 under attribute C 1
Combined with formula (11), the proximity between the attribute values of each scheme of decision-maker E 1 and decision-maker E 2 is: According to formula (12), the average proximity between vectors r 2 1 and r 1 1 is: Similarly,  www.nature.com/scientificreports/ Using formula (13), the proximity between the decision maker E 1 and the other three decision makers on attribute C 1 is: 3. Determine the weight of decision maker E 1 under attribute C 1 In formula (14), take = 0.5 and solve the combination weight ω 1 1 of decision maker E 1 under attribute C 1 .
The combination weights of the other three decision-makers under the attribute C 1 can be obtained by using the similar calculation process. The calculation results are shown in Table 5.
In the same way, the decision-maker weights of all decision-makers on all attributes are calculated, as shown in Table 6.
Step 2: Using the n-dimensional IVIFWM operator in combination with the determined weight of the decision-maker, the individual decision matrix of the decision-maker is aggregated according to formula (19) to get a group decision matrix R = (r ij ) m×n , as shown in Table 7.
Step 3: According to formula (17), the weight w j of each attribute is calculated as follows: The relative weight w ′ j of attributes calculated by formula (18) and Step 4 is: Step 4: According to the dominance matrix under attribute C j condition and referring to the value range of risk aversion coefficient, θ = (0.3, 0.8, 1.0, 2.5, 3.0, 4.0).   Table 6. Decision makers' weights on all attributes.  Table 7. Interval-valued intuitionistic fuzzy group decision matrix. www.nature.com/scientificreports/ For the attribute C 1 , when θ = 1 , the dominance matrix ϕ 1 (A i , A k ) of the scheme A i relative to the scheme A k under the attribute C 1 is calculated by the formula (20) according to the D HH distance measure formula.
In the same way, the dominance matrices of the other attributes ( C 2 ,C 3 ,C 4 ,C 5 ) can be obtained, which is not described in detail in the operation process due to limited space.
Step 5: According to formula (21), the overall superiority of the scheme A i relative to the scheme A k is obtained.
Step 6: According to formula (22), the comprehensive evaluation value of each scheme is obtained, and the schemes are sorted.
It is concluded from Table 8 that when the risk aversion coefficient θ takes different values, the ranking of the schemes is A 2 ≻ A 4 ≻ A 1 ≻ A 3 . According to the ranking of the alternative schemes under different psychological indicators or risk preference degrees, and in combination with the stability presented in the decision-making process, the stable optimal scheme can be determined as A 2 .
As the difference of decision makers' risk preference in practical decision-making problems is objective and cannot be ignored, it is very reasonable and necessary to introduce risk aversion coefficient θ to quantify the degree of risk preference of decision makers in this study. On the premise of the determined risk aversion coefficient θ , selecting any interval-valued intuitionistic fuzzy distance measure formula based on Hausdorff measure in the calculation does not affect the ranking of the alternatives, the selection of the distance measure formula has weak correlation with the determination of the optimal scheme, and the decision result has certain robustness.
If the D HE distance measure algorithm is adopted, the ranking of alternatives will also change obviously with the increase of the θ value, that is, the change of decision makers' attitude towards risks, the risk attitude of decision makers directly affects the whole decision-making behavior. When θ value is small, alternative A 3 is better than A 1 ; With the increase of θ , alternative A 1 is superior to A 3 . In the actual selection decision, if the supplier has unstable economic conditions and institutional sustainability 43 and weak risk resistance, a smaller θ value can be selected in the calculation using this model. On the contrary, a larger parameter θ can be selected. This study shows that the emergency decision-making results under different risk attitudes can be obtained by changing the risk aversion coefficient θ to make the model more applicable. The risk attitude of decision makers is taken into account when sorting the schemes, and different decision options are provided for decision makers with different risk attitudes to better meet people's decision-making needs and enhance the flexibility of decision-making methods. Under the premise of maximizing the profits of supply chain enterprises, the research content herein provides ideas for formulating more accurate cooperation contracts to achieve win-win cooperation among supply chain members. As the research is based on the complexity of human decision-making behavior and risk attitudes, the drawn conclusions will be closer to the reality.
The actual application of a construction machinery manufacturing company in the central and southern China region shows that the optimal supplier well matches the actual supplier selection results. It is an enterprise with a good operation state in the construction industry and a long-term and stable cooperative relationship has  www.nature.com/scientificreports/ been established with the machinery manufacturing company. This shows that this model is feasible to apply in actual supplier selection, which can provide reference for decision makers in supplier selection. Thus it is clear that the method proposed in this study the related research results 42,44,45 has the following advantages over: 1. The method proposed in this paper is compared with the group decision-making method based on simulation and IVIFS proposed in reference 42 , and the ranking result obtained by the combination of simulation and IVIFS is A 2 ≻ A 3 ≻ A 4 ≻ A 1 . It is obvious that the optimal scheme obtained by the method proposed in this study is the same as that obtained by the group decision-making method in reference 42 , that is,A 2 , but the order of each scheme is still different, mainly because the risk preferences of different decision-makers are taken into account in the TODIM method used in this paper, the influence of hesitation on decision-making results is quantified from the perspective of decision-makers, so the decision-making result is more in line with the actual decision-making situation. 2. It is more reasonable and practical to use the method proposed in this paper to determine the weight of decision-makers than the extended TOPSIS method used in reference 44 , because the weights of decisionmakers determined by the extended TOPSIS method are the same for different attributes, and those by the method proposed in this paper are different for different attributes. In actual decision-making, since different attributes may involve different fields, and each decision maker is obviously only good at some fields but not all fields, it is more reasonable that decision makers have different weights under different attributes. In addition, the existing research results use similarity or proximity to determine the weight of decision makers, but the similarity is not necessarily high under great proximity. Therefore, this paper comprehensively uses similarity and proximity in the process of determining the weight of decision makers to better make full use of decision information, so that the determined decision maker weights are closer to the actual weights. 3. Compared with the extended TOPSIS method in literature 44 , the method proposed in this paper has strong distinguishability, which can help us determine the unique optimal solution for different risk preferences of decision makers. Nevertheless, the extended TOPSIS method may result in close comprehensive evaluation value of each scheme, which weakens the distinguishability of this method and makes it difficult to distinguish the schemes. 4. Compared with the method in reference 45 , the focus points of the two methods are obviously different.
The method proposed by the former constructs the interval number fuzzy preference relation based on the unknown function and the fuzzy preference relation, and is applied to the MADM problem with only one decision maker, while the method proposed in this paper is applicable to the MAGDM problem with IVIFN whose attribute weight information is completely unknown. In addition, in the former method, the decision maker only needs to evaluate the preference relationship between schemes after pairwise comparison, while in the present method, each scheme is evaluated with respect to each attribute.
The paper improves the determination method based on decision maker's weight, and determines the decision maker's weight based on the fuzzy degree and uncertainty degree of the judgment information provided by the decision maker, which largely avoids subjectivity and authoritative monopoly. At the same time, considering the influence of different decision makers' risk preference degrees on the decision-making results, to avoid excessive weighting in an interval-valued intuitionistic fuzzy matrix, a measurement coefficient is introduced for the decision makers' risk preference degree. Decision makers can choose the corresponding decision-making results according to their own risk preferences, which can add more choices for decision-makers and increase the flexibility of the decision-making process. To sum up, the method proposed herein can not only provide the ranking results of decision makers in different risk preference styles, but also provide comprehensive decisionmaking results. This method has the advantages of objective and comprehensive consideration of problems and strong operability.

Conclusions
With the rise of outsourcing and procurement globalization, supplier selection has become the focus of the current enterprise supply chain management because high-quality suppliers can effectively reduce the procurement cost of enterprises, gradually increase the flexibility of supply chain management and continuously improve the core competitiveness of enterprises. Under the environment of global procurement, outsourcing of non-core competence, internet and e-commerce, the decision-making of supplier selection is becoming more and more important and complex, and has always been the attention of the theoretical circle and the actual management workers. In this study, by considering the similarity and proximity at the same time, i.e. considering the similarity between the individual decision matrix and the group decision matrix of the decision maker, and referring to the proximity between the decision matrix of a single decision maker and the comprehensive decision matrix composed of other decision makers, a new method is proposed to determine the decision maker's weight for each attribute, which is more objective, comprehensive and persuasive in the discussion of the problem. Finally, an example of mechanical equipment supplier selection decision-making shows the feasibility and effectiveness of the selection evaluation criteria and MAGDM method, which can provide reference for decision-makers to select suppliers, and also provide an effective scientific method to solve the problem of intuitionistic fuzzy MAGDM. In the practical application process, it is necessary to establish relevant databases sand accumulate more experience. Appropriate model parameters can be selected to improve the performance of the group decision-making model of the IVIFS and TODIM method. The method proposed herein can more objectively and reasonably solve the interval-valued intuitionistic fuzzy MAGDM problem when the attribute weight and the decision maker's weight information are completely unknown. It can meet the needs of decision makers with different risk preferences, www.nature.com/scientificreports/ enrich and develop fuzzy group decision-making theory and method, broaden the application scope of group decision theory, thus providing decision makers with more practical and feasible decision-making reference.

Data availability
The data used in this paper are from the raw data of references 42,44 , this paper uses public data (unclassified data) for example analysis, in order to verify the feasibility and effectiveness of the proposed model and method.